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Selective harmonic elimination (SHE) is an efficient method to eliminate 
low-order selected harmonics. However, due to nonlinearity in the problems, 
many optimization techniques give unsatisfied performances in finding 
optimum switching angles for the SHE. This paper proposes a modified 
moth-flame optimization algorithm to eliminate selective harmonics in 
cascaded multilevel inverters. The optimization algorithm is employed to 
find sets of optimum switching angles for cascaded 5-level, 7-level, and 9- 
level inverters. The results have shown that modified moth-flame 
optimization is beneficial in finding optimum switching angles. It performs 
better than moth-flame optimization (MFO) and differential evolution (DE) 
algorithms. The optimum switching angles are applied to generate switching 
pulses for a cascaded 9-level inverter to demonstrate the algorithm’s 
accuracy. As a result, the low-order harmonics are entirely removed from the 
ac output voltage of the inverter. 
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1. INTRODUCTION 

Multilevel inverters have attracted the attention of many researchers due to the need for high power 
applications in industries [1]. In high power applications, the switching frequency of inverters is limited by 
switching losses and electromagnetic interferences. As the switching frequency and switching losses are 
decreased, the efficiency of inverters is increased significantly [2]. 

The advantages of the multilevel inverters are a) they can produce almost sinusoidal waveform 
output voltage with very low total harmonic distortion (THD), b) they can operate at low switching frequency 
[3], and c) they are applicable for high power and high voltage systems. There are commonly three types of 
multilevel inverter topologies: flying capacitor [4], diode-clamp [5, 6], and cascaded multilevel inverter [7]. 
Among the topologies, the cascaded multilevel inverter becomes more attractive for medium and high 
voltage renewable energy systems due to its modularity and simplicity [8]. Another advantage of the 
cascaded multilevel inverter is that it requires less number of components than diode-clamp and flying 
capacitor inverters. 

One of the crucial issues in multilevel inverters is to eliminate harmonics of the ac output voltage. 
Selective harmonic elimination pulse width modulation (SHEPWM) [9-12] is the most popular technique to 
improve the output voltage waveform of multilevel inverters. The SHEPWM technique can be utilized to 
reduce both the switching frequency and THD value of the multilevel inverters. Therefore, it would give a 
significant advantage for high-power multilevel inverters to operate with low switching frequencies [13]. 
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However, accurate estimating of switching patterns plays a crucial role in SHEPWM inverters. The 
main problem associated with the SHEPWM technique is the solution of nonlinear equations under various 
operating conditions of the inverters [14, 15]. 

There are two main methods for solving a set of SHEPWM transcendental equations: deterministic 
[16] and stochastic [17, 18] approaches. Optimization techniques, especially stochastic metaheuristic 
algorithms [19], have become the exciting tools to deal with complex optimization problems in SHEPWM 
inverters. However, due to nonlinearity and multiple variables existed in the issues, many optimization 
algorithms give unsatisfied performances for their premature or slow convergence. 

In this paper, a modified MFO is employed to find sets of optimum switching angles for SHEPWM 
cascaded 5-level, 7-level, and 9-level inverters. The optimization problem for the cascaded multilevel 
inverters is computed by using Matlab. The pulses generated from optimum switching patterns and power 
circuits are built by using Cadence-Pspice. Finally, results are obtained and verified in simulations for a 
cascaded 9-level inverter. 


2. CASCADED MULTILEVEL INVERTER 

A cascaded multilevel inverter consists of N single-phase H-bridge inverters with separate N dc 
sources. The single-phase H-bridge inverters are connected in series, as shown in Figure 1. It also shows that 
the cascaded multilevel inverter has equal dc sources, 


Vacı = Vace = + Vaen = Vac (1) 


Each single-phase H-bridge inverter requires one dc source, or in other words, the cascaded (2N+1)- 
level inverter requires N dc sources [20, 21]. However, the cascaded multilevel inverter can also be supplied 
by a single DC source [22, 23]. The multilevel inverter could produce medium- to high-voltage output from 
low-voltage input. The ac output voltage is synthesized by N dc voltage sources connected to the individual 


H-bridge inverters. 
ic S13 
Vaci 
| Sis S2 
1 S23 
































S2 
Veo = 
Vac3 == 
Ven 
Sna Sno 








Figure 1. The power circuit of cascaded (2N+1)-level inverter 
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The summation of voltages generated by each H-bridge inverter produces a staircase output voltage 
waveform (as seen in Figure 2). Each single-phase H-bridge inverter generates +Vac, 0, and —Vac output, and 
the resulting output voltage of a cascaded (2N+1)-level inverter ranges from —NVac to NVac. 

Figure 2 presents the output voltage waveform of the cascaded (2N+1)-level inverter, which can be 
expressed in Fourier series as follows: 


N 
v(at)= DV, sin(nat) (2) 
n=1,3,5.... 
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Figure 2. The output voltage waveform of a cascaded (2N+1)-level inverter 


The magnitude of harmonic components (including the fundamental) for equal dc sources is given by 


Vael E 
V, = Mal Sy, coin) (3) 
AM | k=l 
The main goal of the SHEPWM is to eliminate (N-1) harmonics from the waveform. The remaining 
equation is used as the magnitude of the fundamental component at the desired value, V;. Based on (3), the 
simultaneous equations for determining N switching angles of the SHEPWM is given by 


ays! 5 sos ie 
fiC ay Pee pa) = El 
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= cos30,)|=e 
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f:(0)= V, = 57M | È e0600 |= es 
Von-ı 4 N 
fonn (0) = a a = | È costa -Dd |=sen-n (4) 


where fs, fs, ..../(2w-1) is the normalized magnitude of harmonics (respect to the fundamental), and M = V/Va- 
is the modulation index. The objective function of the SHEPWM problem is defined as 


E= leil + lezl + les] ++ [Econ-ay| (5) 
The optimum switching angles are obtained by minimizing (5) and must satisfy the following constraint: 
0° < 0, < 0, < ++ < OQn-1) < 90° (6) 


The first of (4) guarantees the desired fundamental component, and the others are utilized to ensure 
the elimination of 3", 5", 7%, 9", ..., and (2N-1) ® harmonics. 
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3. MODIFIED MOTH-FLAME OPTIMIZATION ALGORITHM 

The moth flame optimization (MFO) algorithm uses moths as search agents that move around the 
search space [24]. Flames are considered as the best position of moths in the search space. Each moth 
searches around a flame and updates the discovery as its best solution. 

In solving the SHE problems, the optimum switching angles are the best position of moths in the 
MFO algorithm. The initial parameters of the MFO are the moth population N, the number of variables D, the 
lower bound l, the upper bound u», and the maximum iteration Max. The lower and upper bounds of 
variables are defined as follows 


lb; = [lb, lb) .. lbp], ubj =[ub; ub, .. ubp] (7) 
The moth population is randomly generated by 

M,j;"" = lbj + rand;.(ub; — lbj), i = 1,2,...,N; j =1,2,...,D (8) 
The moth population of the MFO can be presented as 


My, Miz. Mip 
Mz, M22 Mop 


M Ove M31 M32 a M3 (9) 


ij 
My Mya Man 

The fitness value of the moths is carried out by using 

MES? = $M) i= 1,2., N; j = 1,2., D (10) 

The fitness value can be presented as 


MF, 
MF, 
MF;® = | MF; (11) 


MFy 
During the iteration process, the flame fitness (at Z = 1) is sorted of initial fitness values of the 


moths, whereas the flames are sorted according to their fitness values. The flame fitness is sorted from the 
best to the worst values, 


FF, 
FF, 
FF,” = sorted(MF,) = | FF; (12) 


FFy 
Fa View Ke 


Poy F22% Fob 
Fj =| Fai Faz Fap (13) 


fig, Bae Fes 
Updating the moth population is using a logarithmic spiral model as follows 
M; jl? = $,;,°.e©. cos(2nt) + Fj? (14) 
where S is the distance between the moths and the flames, b is a shape constant of the spiral path (-1 < b < 1), 


t is a control parameter to maintain the distance between moth and flame in the spiral path (r < t < 1), and ris 
a convergence constant which is decreasing from -1 to -2 to estimate the values of t. 
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The distance between the moths and the flames is defined as 
Sij = abs(F;j® — M;;j®), i= 1,2., N; j = 1,2, n, D (15) 
whereas r and t are as follows 
r= —(1+—~),t=rand.(r-1) +1 (16) 


The position of updating moths could eventually degrade after passing iterations. In order to ensure 
that the moths will meet the flames, the moths should update their positions with respect to the best flames at 
the final step of the iteration. A control parameter necessary for this mechanism is 





Pa _ (N-19) 
FN = round (N — ) (17) 
According to (17), the moth population is updated by using 
M; jP = $,;°.e©9. cos(2nt) + Fj if I < FN 
MU = $,;,©.e©9. cos(2nt) + Fry, j© if 1 > FN (18) 
while updating the flames: 
MF; (I+1) 
BR = Sorted | $ p hi= toa N (19) 
FF, 
The flame corresponding to its fitness value is 
Fa Fe Fy 
aa F22% Toy 
ie! = F34 F 32 any Fz p (20) 
Fyi Fuze Fup 
Finally, the best position of the moths and their fitness values are selected in here as 
EFV tD = FARO 
Ft P SR ol S G DeD (21) 


Updating the moths and the flames processes are repeated until FFbest(I+1) meets the criterion € as 
defined in (5) and/or I = Max. 


4. RESULTS AND ANALYSIS 
4.1. Switching patterns 

Achieving the fundamental component at a desired of M and suppressing selected harmonics to zero 
as describing in (4) are very difficult to be solved numerically. In this work, a modified MFO algorithm is 
applied to alleviate the computational problems. The optimum switching angles for cascaded 5-level, 7-level, 
and 9-level inverter are explored by using the modified MFO for all values of M. As the results, the optimum 
switching angles for the cascaded 5-level, 7-level, and 9-level inverter are presented in Figure 3, 4, and 5, 
respectively. All optimum switching angles must meet the requirements of (6), as mentioned earlier. Thus, 
there are no possible optimum switching angles within M < 1.1027 for a cascaded 5-level inverter. Using the 
same constraint which must satisfy (6), there are also no optimum switching angles within M < 2.0974 for a 
cascaded 7-level inverter and within M < 3.0930 for a cascaded 9-level inverter. It is noted that not all M has 
its solutions for the optimum switching angles since the optimum switching angles will only be applicable if 
they fulfill the requirement as defined by (6). 
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Figure 4. Optimum switching patterns for cascaded 7-level inverter 


Figure 5. Optimum switching patterns for cascaded 9-level inverter 
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4.2. Modified MFO performance 

The modified MFO algorithm is compared to MFO and DE algorithms in order to verify its 
performance. These algorithms are analyzed by using the same objective function, as mentioned in (4-6). The 
approach using is to stop iterations when the searching processes have reached £ < 107. Twenty independent 
experiments are carried out for each algorithm. As the results, the average number of iterations required by 
the modified MFO is 115 iterations, while the MFO and the DE need 131 and 182 iterations, respectively. 
Comparison to both MFO and DE algorithms, it is revealed that the modified MFO gives better performance 
in finding optimum switching angles (as seen in Figure 6). The modified MFO is convergence faster than 
MFO and DE algorithms. 
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Figure 6. Fitness value curves for ¢ < 107 


4.3. Cascade 9-level inverter 

In order to generate output ac voltage from a cascaded multilevel inverter, the parameters used in 
simulations are: cascaded 9-level, Va = 100V, M = 3.2, and f = 50 Hz (T = 20 ms). The set of optimum 
switching angles found for M = 3.2 are: 6; = 10.8169°, 02 = 26.3546°, 03 = 53.0106°, and 64 = 88.0910°. The 
time delay and pulse width during the PWM control of the inverter are derived from the optimum switching 
angles. The switching patterns are presented in Figures 7 and 8. 

The simulation results at a modulation index M = 3.2 are given in Figure 9. From the frequency 
spectra, it can be seen that the magnitudes of low-order harmonics such as the3", 5, and 7 have been 
omitted from the waveform entirely. 

The residual higher-order harmonics are: Vo/V; = 3.6665%, Vii/V1 = 4.4569%, Vi3/Vi = 4.4684%, 
Vis/Vi = 0.9545%, Vi7/Vi = 3.3251%, Vio/Vi = 4.0786%, V21/V; = 0.3430%. The measured value of V; is 
319.843 V. It matches with the calculated V; for M = 0.8, which is V; = 310 V. Each higher-order harmonic 
component is lesser than 5% with total harmonic distortion (THD) of 9.95%. The THD of output voltage can 
be reduced by applying an LC passive filter to minimize the high-order harmonics [25]. 
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Figure 7. Switching pulses for a cascaded 9-level inverter: S11, S12, S21, S22, S31, S32, Sai, S42 
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Figure 8. Switching pulses for a cascaded 9-level inverter: S13, S14, S23, S24, S33, S34, S43, S44 
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Figure 9. The output voltage waveform and the corresponding FFT analysis 


5. CONCLUSION 

A modified MFO algorithm has been proposed in order to overcome computational difficulties of 
SHEPWM cascaded multilevel inverters. The modified MFO algorithm provides precise computation of 
optimum switching angles for cascaded multilevel inverters. It is faster convergence than both MFO and DE 
algorithms. Simulation results are presented for a cascaded 9-level inverter in single-phase configuration. 
The results also show that the low-order harmonics are utterly eliminated from the ac output voltage 
waveform. It reveals that the modified MFO is an effective algorithm to solve SHE problems of cascaded 
multilevel inverters. 
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